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In a model with a gauge group Gsm ® U{l)x, where Gsm = 517(3)0 ® SU{2)l ® U{1)y is 
the standard model gauge group and U{l)x is a horizontal local gauge symmetry, we propose 
a radiative generation of the spectrum of quark masses and mixing angles. The assignment of 
horizontal charges is such that at tree level only the third family is massive. Using these tree 
level masses and introducing exotic scalars, the light families of quarks acquire hierarchical masses 
through radiative corrections. The rank three quark mass matrices obtained are written in terms of 
a minimal set of free parameters of the model, whose values are estimated performing a numerical 
fit. The resulting quark masses and CKM mixing angles turn out to be in good agreement with the 
experimental values. 
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I. INTRODUCTION 

Although for several years a great effort has been done to shed some light on the mystery of the fermion masses it 
still is one of the outstanding puzzles of particle physics. There have been different approaches to explain the mass 
hierarchy, the fermion mixing and their possible relation to new physics. A good review covering widely these topics 
has been presented in Q . In this work we will restrict our study to that of a new horizontal symmetry and the derived 
radiative corrections. 

A possible answer to why the masses of the light quarks are so small compared with the electroweak scale is that 
they arise from radiative corrections Q , while the mass of the top quark and possibly those of the bottom quark and 
of the tau lepton are generated at tree level. This may be understood as a consequence of the breaking of a symmetry 
among families (a horizontal symmetry). This symmetry may be discrete or continuous!^. Here we consider the 
case of a continuous local horizontal symmetry, a U{l)x gauge group broken spontaneously. We limit our calculation 
of masses to those of the quark sector, insisting that at tree level only the top and bottom quarks acquire mass. 
Instead of assuming a texture for the quark masses from the beginning, we carry out a one loop and a two loop 
calculation of the mass matrices in terms of some parameters which are the tree level top and bottom quark masses, 
one Yukawa coupling and the entries in the mass matrices of the scalar bosons that participate in the loop diagrams 
which contribute to the quark masses. 

This paper is organized in the following way: In Section II wc describe explicitly the model. Section III contains the 
analytical calculations, while Section IV is devoted to a numerical fit of our equations. The conclusions are presented 
in section V. 



II. THE MODEL 



We assume only three families, the Standard Model (SM) families, and we do not introduce exotic fermions to 
cancel anomalies. The fermions are classified as in the SM in five sectors / — q,u,d,l and e, where q and / are the 
SU(2)l quark and lepton doublets respectively and u,d and e are the singlets, in an obvious notation. In order to 
reduce the number of parameters and to make the model free of anomalies, we demand that the values X of the 
horizont al charge satisfy the traceless condition [|| 
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X(/,) = 0,±5/, (1) 
where i = 1, 2, 3 is a family index, with the constraint 

51-251 + 51 = 5^-51 (2) 

Eq. (nl) guarantees the cancellation of the [U(1)h]'^ anomaly as well as those which are linear in the \]{1)h hypercharge 
{[Su\z)cYU{1)h, [SU[2)LfU{l)H, [Gmt;] 2 [/(I) ^ and [U{1)yYU[1)h)- Eq. (||) is the condition for the cancellation 
of the U(l)y [U(l)ij]2 anomaly. A solution of Eq. which guarantees that only the top and bottom quarks get 
masses at tree level is given by ( "doublets independent of singlets" , see Ref. |^ ) 

5i=5q = ±^^ 5., = 5d = 5, = ±5. (3) 

To avoid tree level flavor changing neutral currents, we do not allow mixing between the standard model Z boson and 
its horizontal counterpart. Consequently the SM Higgs scalar should have zero horizontal charge. As a consequence, 
and since we insist in having a non-zero tree-level mass for the top and bottom quarks, the horizontal charges of these 
quarks should satisfy 



- X(<73) + X{u:,) = 0, -X(g3) + X{dz) = (4) 

in order for the Yukawa couplings in Eq. (|^) to be invariant, but then Eqs. (|l|) and (|^) demand that they vanish, 

X(u^) = Xiqs) = Xids) = 0, (5) 

which in turn implies ^(^3) = ^'(63) = (this defines the third family). The assignment of horizontal charges to the 
fermions is then as given in Table 1. The SU{3)c ® SU{2)l ® U{1)y quantum numbers of the fermions are the same 
as in the Standard Model. 



Scctor\ Family 
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±5 
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TABLE I: Horizontal charges of fermions. 
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TABLE IL Quantum numbers for scalar fields, C denotes the dimension of the representation under the SU{3)c color group. 



To generate the first and second family quark masses radiatively we must introduce new irreducible representations 
(irreps) of scalar fields, since the gauge bosons of G = Gsm ^ U{l)x do not perform transitions between different 
families. Families are of course distinguishable (non degenerated) only below the scale of the SM symmetry breaking, 
when they become massive. 



3 



Looking for scalars which make possible the generation of fcrmion masses in a hierarchical manner, we divide the 
irreps of scalar fields into two classes. Class I (II) contains scalar fields which get (do not get) vacuum expectation 
value (VEV). 

A proper choice of scalars should be such that no VEVs are induced, through couplings in the potential, for scalars 
in class II. In the model considered below scalars of class II have no electrically neutral components, so they never 
get out of its class. In our model we introduce two irreps of scalars of class I and six irreps of scalars of class II, 
with the quantum numbers specified in Table Notice that we introduce just the minimum number of scalars of 
class I; that is, only one Higgs doublet of weak isospin to achieve the Spontaneous Symmetry Breaking (SSB) of the 
electroweak group down to the electromagnetic U{1)q, and one 5*1/(2)^ singlet 02 used to break U{l)x- In this way 
the horizontal interactions affect the p parameter only at higher orders. 

With the above quantum numbers the quark Yukawa couplings that can be written may be divided into two classes, 
those of the D type which are defined by Fig la, and those of the M type which are defined in Fig lb. The Yukawa 
couplings can thus be written as Ly = Ly^ + -^Ym j where the D Yukawa couplings are 



Lyo = Y^qLs4>lU]i^ + F'^gLa'/'l'^-Ra + 'i-c-, (6) 

with (j) = *f20*, while the M couplings compatible with the symmetries of the model are 

+U2RCcj)jUiii + u'^j^C(j)7U3R + u'^RC(j)sU2R] + h.c. (7) 

In these couplings C represents the charge conjugation matrix and a and /3 are weak isospin indices. Color indices 
have not been written explicitly. By simplicity and economy we have assumed only one Yukawa constant Yj for all 
the M couplings. Notice that (p^afS} is represented as 



_ / 0-4/3 0-1/3 \ 
"P^^y 0-1/3 02/3 j 

where the superscript denotes the electric charge of the field. The same applies for 04. 



(8) 



Scalar fields which are not SU(2)i doublets do not participate in D type Yukawa terms, they however contribute 
to the mass matrix of the scalar sector and in turn determine the magnitude of the radiatively generated masses of 
fermions, as we shall see below. 

The most general scalar potential of dimension < 4 that can be written is 



- y(0o ^Y.p^\cj>,\^ +Y, K I -^^ n <t>3 ? +v^A<t>i<t>^<i>i + ^~3i0i 0^0301 

+r?410l04<^401 + Vll't'l^ 4>iMl + X! '^'^ I '^i'^J' 1^ +iPl4>l4>6(t>2 + 

*,i#l,2 

P207'/'802 + Ai0|^0i 03{a/3}0^* + A2 0701 " 03{q/3} 0^ + 
A3Tr(0|^04)02 + A40506'?!'702 + A5050|j0|08 + h.c), (9) 

where Tr means trace and in | 0^ p= 0j0i an appropriate contraction of the SU{2)l and S't/(3)c indices is 
understood. The gauge invariance of this potential requires the relation A = 25. 

Now we proceed to describe the mechanism that produces the quark masses. In general we could have contributions 
of two types as depicted in Fig ^ In the present model however, we have only the diagrams of Fig ^ for the charge 
— 1/3 quark mass matrix elements and similar ones for the charge 2/3 sector (these type of diagrams were first 
introduced in ) ; in those diagrams of Fig |^ the cross means tree level mixing and the black circle means one loop 
mixing. The diagrams in Fig ^ y|^ should be added to the matrix elements (1,3) and (3,1), respectively. In the one 
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fi"- fsL f3R fjR 



FIG. 1: Generic diagrams that could contribute to the mass of the light families, (a) D type couplings are represented with 
vertices where one fermion is incoming and the other one is outgoing, (b) the M type couplings are represented with vertices 
where both fermions are incoming or outgoing. 



d IL d 2L d 2R d IR d 2L d 3L d 2R d IR d 3L d 3L d 2R d IR 
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FIG. 2: Mass matrix elements for d quarks 



loop contribution to the mass matrices for the different quark sectors only the third family of quarks appears in the 
internal lines. This generates a rank 2 matrix, which once diagonalized gives the physical states at this approximation. 
Then using these mass eigenstates we compute the next order contribution, obtaining a matrix of rank 3. After the 
diagonalization of this matrix we get the mass eigenvalues and eigenstates (A quark mass mechanism with some 
similar features to the one proposed here is given in 



d 3L d 2L d 2R d IR d IL d 2L d 2R d 3R 

(a) (h) 



FIG. 3: Extra diagrams that contributes to mass matrix elements (a)(l,3) and (b) (3,1) 
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The VEVs of the class I scalar fields are 



< 02 >= V2, 



(10) 



and they achieve the breaking 



GsM ® Uil)x GsM SU{3)c ® t/(l)c 



(11) 



The scalar field mixing arises after SSB from the terms in the potential that couple two different class II fields to 
one of class I. After SSB the mass matrices for the scalar fields of charge 2/3 (04,03,05,06) and 4/3 (04,03,07,08) are 
written, respectively, as 



^2/3 = 
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(12) 



P2V2 H 



4 I 



where from Eq. (|9|) if = wf = ^| + \i\v\ + Ai2w| and ^ rjuvj. 

Notice that due to the scalar mixing in all the loop diagrams of Fig || and |[ the divergences in each one of these 
diagrams cancel as is physically expected, giving rise to finite contributions to the quark mass matrices. 

Explicitly, the non vanishing contributions from the diagrams of Fig ^ to the mass terms dmdjL^l^^ + h.c. read at 
one loop 



.(1) 
■^22 



3™ 



(0) 



Y, 



(0)^ 



(13) 



4y = ^ C/2feC/4./(M,, m(°)), 

k 



(0) Y,^ 



(14) 



(15) 



where w^''^ is the tree level contribution to the b quark mass, the 3 is a color factor, U is the orthogonal matrix which 
diagonalizes the mass matrix of the charge 2/3 scalars, 

(04,03,05,06)"^ = [/(cri,cr2,(T3,0-4)'^, 

where ai are the eigenfields with eigenvalues Mi, and 

1 



/(a, 6)^ 



which is just a logarithmic contribution when ^ b^. The resulting second rank mass matrix at this level is thus 
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^32 


(0) 

ml 



(16) 



At effective two loops we obtain the following expressions: 
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k.i 



^12^ - 3^^mW(F«)3.(T/ij,))2,C/ifeC/3fe/(Mfe,m|^)), (18) 

k.i 



43^=3^^mf)(V^«)3.(Fij,^)2.C/2fct/3fc/(M,,mf)) (19) 



3^ E "^rH^il^)2.(V;(H^)2.C/uC/3fe/(Mfe, mf^ 



IGtt 

k,i 

S^i^ = 3^5^mW(T/«)2.(l^ijj^)3,C/2feC/4fe/(Mfe,m«), (20) 

k^i 

4f =3^E"^^^'^(^ii)2^(^i«)3»t^2fe^3./(M,,m«) (21) 
+3l£E"^r^(^iL^)2.(V;,^H^)2.C/2.C/4fe/(Mfe,mW), 

where the k (i) index goes from 1 to 4 (from 2 to 3), V^^ and vj^ are the unitary matrices which diagonahze M^^'' 
of equation ([l6|) and m|^'' are the eigenvalues. Therefore at two loops the mass matrix for d quarks becomes: 



Mf = 



/ v(2) y(2) y(2) 
^11 ^12 ^13 

E^'i' 1 . (22) 



For the up sector the procedure to obtain the masses is completely analogous. That is, the mass terms for the up 
sector come from graphs like those in Fig. ^ and but replacing the 04, (/'3, (f)^ and scalar fields by (f>4,4>3, (l>j and 
(/)8 and the quarks di by the quarks Ui. 

The CKM matrix takes the form 

V^CKM = (K!lVW)tFi^V«, (23) 

where the unitary matrices V^j^ and V^j^ diagonalize mI^'\ and V^j^ and V^^l diagonalize with an analogous 

notation used for the down sector. 

It is important to mention here that the textures, particularly the zeros in the scalar and quark mass matrices (Eqs. 
[l^ and [2^ ) are not accidental neither imposed; they are just a direct consequence of the mass mechanism that we are 
introducing and of the gauge symmetry of the model. 
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III. NUMERICAL EVALUATION 



A. Experimental values 

Since quarks are confined inside hadrons, their masses can not be directly measured. So, the quark mass parameters 
in the SM Lagrangian depend both on the renormahzation scheme adopted to define the theory and on the scale 
parameter /i where the theory is being tested. In the limit where all quark masses are zero, the SM has an SU (3) l ® 
SU{3)ii chiral symmetry which is broken at an scale ~ IGeV. To determine the quark mass values one must use 
SM perturbation theory at an energy scale /i >> where non perturbative effects are negligible. 

For illustration, the allowed ranges of quark masses^ in the modified minimal subtraction scheme (MS) are|0: 

muil.GeV) = 2-6.8MeV. 
nidil-GeV) 4- 12 MeV. 
msil.GeV) = 81- 230 Afey. 

mc(TOc) = 1-1 - 1-4 GeV. 

mb{mb) = A.l-AAGeV. 
mt{Exp) = m.%±b.2GeV 

To get the relative magnitude of different quark masses in a meaningful way, one has to describe all quark masses 
in the same scheme and at the same scale. In our analysis we are calculating the quark masses at an energy scale 
/i„i such that Mz < ^J■m < Mx — V2, where Mx is the mass scale where U{l)x is spontaneously broken. Since in 
our model there is no mixing between the Standard Model Z boson and its horizontal counterpart, we can have V2 
as low as the electroweak breaking scale. For simplicity, let us assume that our calculations are meaningful at the 
electroweak breaking scale and from the former values for the quark masses let us calculate, in the MS scheme, the 
quark masses at the rrit scale |^ and at the Mz scale Those values calculated in the references cited, are presented 
in tables III and |^ respectively. 



On the other hand, the CKM matrix elements are not ill defined and they can be directly measured from the 
charged weak current in the SM. For simplicity we assume that they are real, and as discussed in Ref B, they are 



almost constant in the interval Mz < /i < a few TeV. Their current experimental value p] are given in the Tables [II 



and IV 



B. Evaluation of the parameters 

In order to test the model using the least possible number of free parameters, let us write the scalar mass matrices 
in the following form: 
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a_ 
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d- 





(24) 



Using the central value of the CKM elements in the PDG book|^ and the central values of the six 
quark masses at the top mass scale P, we build the function in the ten parameter space defined by 



(a+, a_, 6, c+, c_, d+, d_, Yf , mjj"'', mj"''), where m}^' and m^' are the tree level quark masses for the bottom and 
top quarks respectively. Expressions for the eigenvectors and eigenvalues of the mass matrices involved in the nu- 
merical evaluation were obtained using MATHEMATICA, and the function was minimized using MINUIT from 
the CERNLIB packages^; both Monte Carlo and standard routines were used in the minimization process. The tree 

level masses of the top {m^i^) and bottom {m^b^) quarks were restricted to be around the central values ± 10 % in 
order to assure consistency with the assumption that radiative corrections are small. The x^ function presents an 
even symmetry with respect to 5 of the parameters of the matrices in Eq. (24); we find that there are 32 parameter 
domains where the x^ function takes small values. For the extremal points this even symmetry is not an exact one, but 
all the zones have a Yukawa constant of the order of 10 and give masses and CKM matrix elements in good agreement 
with the available experimental values. The numerical results on one of the 32 minima of the ten parameter space 
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RANGE 


INPUT 


BEST FIT 


md{mt) 


3.85-5.07 MeV 


4.46 MeV 


3.63 MeV 


ms(mt) 


76.9-100 MeV 


88.4 MeV 


44.9 MeV 


mb{mt) 


2.74-2.96 GeV 


2.85 GeV 


2.91 GeV 


mu{mt) 


1.8-2.63 MeV. 


2.21 MeV 


2.22 MeV 


mc{mt) 


587-700 MeV 


643 MeV 


841 MeV 


mt{mt) 


159-183 GeV 


171 GeV 


166.5 GeV 


CKMii 


0.9745-0.9760 


0.9752 


0.9761 


CKM12 


0.217-0.224 


0.2200 


0.2179 


CKM13 


0.0018-0.0045 


0.0034 


0.0032 


CKM21 


0.217-0.224 


0.2200 


0.2214 


CKM22 


0.9737-0.9753 


0.9755 


0.9742 


CKM23 


0.036-0.042 


0.0390 


0.0382 




0.004-0.013 


0.0085 


0.0117 


CKM32 


0.035-0.042 


0.0385 


0.0365 


CKM33 


0.9991-0.9994 


0.9992 


0.9992 



TABLE III: Experimentally allowed values for mqijnt) and CKM matrix elements. We show the input and calculated values 
in the context of our model. 





RANGE 


INPUT 


BEST FIT 


md{Mz) 


1.8-5.3 MeV 


3.55 MeV 


3.44 MeV 


ms{Mz) 


35-100 MeV 


67.5 MeV 


39.6 MeV 


mh{Mz) 


2.8-3.0 GeV 


2.9 GeV 


2.9 GeV 


mu{Mz) 


0.9-2.9 MeV 


1.19 MeV 


2.03 MeV 


mc{Mz) 


530-680 MeV 


605 MeV 


793 MeV 


mt{Mz) 


168-180 GeV 


174 GeV 


166.9 GeV 


CKMii 


0.9745-0.9760 


0.9752 


0.9762 


CKM12 


0.217-0.224 


0.2205 


0.2174 


CKM13 


0.0018-0.0045 


0.0036 


0.0030 


CKM21 


0.217-0.224 


0.2205 


0.2215 


CKM22 


0.9737-0.9753 


0.9745 


0.9741 


CKM23 


0.036-0.042 


0.0390 


0.0387 


CKM31 


0.004-0.013 


0.0085 


0.0116 


CKM32 


0.035-0.042 


0.0385 


0.0370 


CKM33 


0.9991-0.9994 


0.9992 


0.9992 



TABLE IV: Experimentally allowed values for m,q{Mz) and CKM matrix elements. We show the input and calculated values 
in the context of our model. 



are shown in Table 0. We use those values (which minimize x^) to calculate, in the context of our model, the fifteen 
predictions for m,q(mt) for q — u,d, c, s, t, b and (CKM)ij for i, j = 1,2, 3. The numerical results are shown in Table 

For the sake of comparison, we repeat the same calculations but now using the central values of the six quark 
masses at the AIz scale B. The numerical results are shown in Table IV 



Let us make two comments: First, the values for the parameters in the scalar field square mass matrices are of order 
lO^^(MeV)^ (see Table ^), so, the scalar physical masses are of order 10^ TeV. Second, the rounding errors allow us 
to take safely up to five significative figures in the masses and in the CKM matrix elements 



As can be seen from Tables III and IV, even under the assumption that the CKM matrix elements are real, the 
numerical values are in good agreement with the allowed experimental results. 
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PARAMETER 


a_ 


h 


c_ 






VALUE 


22.845x1016 


1.850 xlQis 


10.018 xlQis 


-10.227 xlO^s 


13.571 xlQi^ 


PARAMETER 


d+ 


04. 


r 


(0) 


(0) 


VALUE 


13.762 xlQi^ 


150.01x10^6 


13.6 


2.912 GeV 


166.1 GeV. 



TABLE V: Values of the parameters in one of the minima (the values of the scalar mass matrices elements are in (MeV)^ 
units). 



IV. CONCLUSIONS 



By introducing a U{l)x gauge flavor symmetry and enlarging the scalar sector, we have presented a mechanism and 
an explicit model able to generate radiatively the hierarchical spectrum of quarks masses and CKM mixing angles. 
The horizontal charge assignment to particles is such that we do not need to go beyond the known three generations 
of quarks and leptons. Also, at tree level only the t and b quarks get masses. To generate radiatively the masses 
for the light families we have introduced some new exotic scalars. All of these new scalars are charged and color 
non-singlets, so they can not get VEV as is required in the loop graphs. 



Our numerical results are presented in Tables III and IV. Even though we are guessing the U{l)x mass scale, the 



two sets of results do not differ by much and they agree fairly well with the experimental values, meaning that the 
mass scale associated with the horizontal symmetry may be in the range IQOGeV < Mx < l.OTeV^. A closer look 
to our analysis shows that we are translating the quark mass hierarchy to the quotient vxjvi which is the hierarchy 
between the electroweak mass scale and the Horizontal U{X)x mass scale. In this way we demonstrate the viability 
that new physics at the electroweak mass scale, or just above it, may help to explain the long-lasting puzzle of the 
enormous range of quark masses and mixing angles. 

Since quarks carry baryon number B = 1/3, the color sextet scalars we have introduced must have B — —2/3 
(the scalar singlets 4>\ and (\)2 have B=0); in this way Ly^ is not only U(X)x invariant but conserves color and 
baryon number as well. On the other hand, 1^(</') does not conserve baryon number; as a matter of fact, the term 
A4(/)5(/)607'/>2 violates baryon number by two units (Ai? = ±2) and could induce neutron-antineutron oscillations. A 
roughly estimate of such oscillations shows that they are proportional to i'2(M|^M|^M|^)~-'^ ~ Vd~^'^GeV which is 
negligible in principle. Any way, in the worse of the situations, since the offending term does not enter on the mass 
matrix for the Higgs scalars (it is just there), it may be removed in more realistic models by the introduction of a 
discrete symmetry. 

Our results are encouraging; even under the assumption that the CKM matrix is real, and without knowing exactly 
the TJ{X)x mass scale, the numerical predictions are in the ballpark, implying also a value of order 10 for the Yukawa 
coupling Yf, and masses for the exotic scalars being of order 10"^ TeV. Our model presents thus a clear mechanism 
able to explain the mass hierarchy and mixing of the quarks. 

Finally let us mention that in the work presented here, the Higgs scalar used to produce the SSB of the SM gauge 
group down to 5't7(3)c x J7(1)q has zero horizontal charge, and as a consequence the standard Z boson does not mix 
with the horizontal counterpart. However, due to recent interest JTo| on the phenomenology of a Z\ it is worth to 
study the possibility of allowing this mixing in future work. 
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